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ABSTRACT 

The  interactions  between  a  semi-submersible  drone  and 
its  large  towfish,  where  the  tension  in  the  tow  cable  join¬ 
ing  them  is  about  20%  of  the  weight  of  the  drone,  are 
investigated.  Since  standard  cable  dynamics  modelling 
computer  programs  do  not  model  tow- vehicle /tow-body 
interactions,  an  iterative  approach  to  calculating  these 
interactions  is  developed.  Drone  motions  are  modelled 
with  a  non-linear,  six  degree-of-freedom,  underwater  ve¬ 
hicle  simulation  incorporating  pre-defined,  three  compo¬ 
nent,  time  dependent  tow  cable  tensions  applied  at  the 
tow  point.  Cable  dynamics  are  modelled  with  the  non¬ 
linear,  finite  segment,  cable/towflsh  dynamic  simulator 
DYNTOCABS,  which  accepts  pre-defined  tow-point  ac¬ 
celerations  as  time  varying  boundary  conditions.  The  in¬ 
teractions  are  calculated  by  iterating  between  these  two 
programs.  'The  method  is  applied  to  simple  turning  ma¬ 
neuvers  important  to  minehunting  operations. 

KEY  WORD'S:  Cable  Dynamics,  Towfish,  Semi-Submersible. 
Coupled  Motions,  Interactions 

NOMENCLATURE 


D 

diameter  of  turn. 

i 

iteration  index. 

£ 

drone  length. 

m.m1 

drone  mass;  m'  —  m/(|pf3). 

K 

yaw  moment  due  to  yaw  rate. 

K 

yaw  moment  due  to  lateral  velocity. 

N'St 

yaw  moment  due  to  rudder  deflection. 

r,  r' 

yaw  rate:  r'  =  rljU . 

rpm 

revolutions  per  minute. 

t. 

time,  in  a  U-turn,  at  which  the  rudder  is  re-zeroed. 

T' 

-‘v 

lateral  body  reference  frame  tension. 

U 

drone  speed. 

XG,  x'o 

drone  longitudinal  center  of  gravity:  x'a  =  xajl. 

Xt , 

longitudinal  location  of  cable  towpoint  on  drone. 

X,  Y,  Z  inertial  reference  frame  displacement  coordinates. 

Yj  lateral  force  due  to  yaw  rate. 

Yj  lateral  force  due  to  lateral  velocity. 

Yir  lateral  force  due  to  rudder  deflection. 

&r,i  i'h  iteration  rudder  deflection. 

tp  drone  heading. 

4>f  final  drone  heading  after  a  manoeuvre. 

8  drone  pitch, 

i i>  drone  roll. 

Primed  quantities  are  dimensionless.  Forces  are  nondimensioEal- 
ized  with  pU2t2  [2,  and  moments  with  pU2£3/2,  where  p  is  the 
density  of  sea  water. 

1.  INTRODUCTION 

DREA  and  ISER  are  evaluating  issues  facing  the  development  of  a 
Canadian  Remote  Minehunting  System  (CRMS)  for  the  Canadian 
Navy.  The  CRMS  is  an  autonomous,  snorkelling  drone  towing  a 
deployable,  active  towfish  which  houses  a  side  scan  sonar  for  route 
surveying  and  mine  location  on  the  sea  floor.  With  route  survey¬ 
ing,  sonar  images  are  obtained  from  an  area  where  mine  hunting 
is  anticipated,  in  order  to  provide  a  reference  against  which  future 
mine  hunting  images  can  be  compared.  A  high  degree  of  towfish 
stability  is  required  to  get  good  images,  and  the  absolute  location 
of  the  towfish  must  be  known  for  differencing  images  with  those 
from  subsequent  mine  hunting  surveys.  Drone  motions,  whether 
from  waves  or  necessary  maneuvers,  cause  towfish  (sonar)  distur¬ 
bances.  These  disturbances  degrade  the  sonar  image  and  make 
mine  detection  and  classification  difficult. 

The  requirement  for  high  towfish  stability  leads  to  a  need 
for  a  stable  towing  platform.  The  serai-submersible  DOLPHIN 
(Deep  Ocean  Logging  Platform  for  Hydrographic  Instrumenta¬ 
tion  and  Navigation)  Mkl  vehicle  developed  by  ISE  Research  Ltd. 
is  a  proven  stable  remote  platform  for  hydrographic  instrumen¬ 
tation.  It  can  operate  in  up  to  sea  state  5  and  has  successfully 
demonstrated  its  towing  capability  (Preston  and  Shupe.  1993). 
DOLPHIN  Mk2  (Figure  1)  is  currently  a  candidate  drone  for 
the  CRMS.  It  is  8.5  m  long  with  aim  hull  diameter,  has  a 
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Fig.  1  DOLPHIN  Mk2  and  towfish. 


dry  mass  of  4500  kg,  and  uses  a  350  horsepower  engine  to  drive 
contra-rotating  propellers.  The  towfish  considered  in  this  study 
(Figure  1)  is  3.2  m  in  length  with  a  0.5  m  diameter  hull.  It  has 
a  dry  mass  of  445  kg,  an  active  depressor  with  a  2.4  m  span,  and 
four  symmetrical,  independently  active  tailfins. 

A  semi-submersible  minehunting  drone  has  personnel  safety 
and  cost  advantages  over  a  minehunting  surface  ship,  yet  retains 
the  larger  ship’s  stability  in  waves,  its  air  breathing  endurance, 
high  data  rate  real  time  communications  link,  and  accurate  global 
positioning  characteristics.  The  drone  is,  however,  still  subject  to 
surge  motions  in  high  sea  states  and  these  result  in  cyclic  tow  ca¬ 
ble  tension  and,  potentially,  towfish  positioning  variations.  Also, 
the  lower  inertia  of  the  drone  results  in  increased  interactions 
with  the  towfish  and  cable.  The  cable  tensions  with  a  towfish 
depressed  to  100  m  depth  with  300  m  of  scope  can  be  20%  of  the 
weight  of  the  drone. 

Therefore,  there  is  a  need  to  assess  the  impact  of  both  drone 
motions  on  the  tow  system  and  the  tow  system  response  on  the 
maneuverability  of  the  drone.  The  dynamics  of  maneuvering  with 
towed  systems  have  been  studied  by  many  workers  (Chapman, 
1984;  Henry,  1985;  Duvat  and  Large,  1997)  who  have  looked  at 
turning  maneuvers  and  how  to  minimize  their  impact  on  the  ship 
(Le  Guerch,  1987)  or  the  towfish  displacement  (Chapman,  1984). 
Ship  induced  motions  on  the  towfish  have  always  been  of  interest 
(Chapman,  1982)  since  a  ship  in  a  sea  state  inevitably  affects 
the  towfish  tethered  to  it.  Some  work  has  also  been  done  on 
towfish  interactions  with  the  tow  vehicle  (Zhu  and  Li,  1997)  but 


only  at  low  speed.  Cable  dynamics  tools  like  BCABLE  (Wilkie, 
1989)  and  DYNTOCABS  (DYNamics  of  TOwed  CABle  Systems, 
Kamman  et  al.,  1989,  1990,  1996)  are  successfully  used  to  model 
tow  cable  systems.  These  programs  do  not  model  the  tow  vehicle 
because  it  is  generally  much  more  massive  than  the  towfish,  so 
little  interaction  occurs. 

The  emphasis  in  this  paper  is  to  develop  a  simulation  tool  with 
which  to  carry  out  interaction  assessments,  and  do  a  preliminary 
assessment  of  one  required  maneuver.  The  interaction  analysis  is 
best  done  with  a  fully  integrated  drone/tow  system  simulator,  a 
tool  unavailable  to  the  authors.  To  get  preliminary  results  quickly, 
and  assess  the  need  for  a  fully  integrated  simulator,  the  analysis 
was  carried  out  by  iterating  between  two  simulation  packages  that 
were  available:  the  in-house  DREA  Submersible  Simulation  Pro¬ 
gram  (DSSP)  still  in  development,  and  DYNTOCABS.  DSSP  in¬ 
tegrates  the  6  degree-of-freedom  equations  of  motion  (Abkowitz, 
1969)  to  obtain  the  vehicle  state  vector  in  response  to  a  series 
of  commands  and  external  forces.  DYNTOCABS  determines  the 
towfish  and  cable  state  vectors  in  response  to  a  given  towpoint 
trajectory.’ 

Minehunting  operations  consist  of  towing  the  sonar  back  and 
forth  to  scan  consecutive  swaths  of  area.  This  requires  periodi¬ 
cally  reversing  the  direction  of  the  entire  tow  system.  The  U-turn 
is  the  simplest  maneuver  with  which  to  do  this  but  it  results  in 
large  towfish  positional  displacements  and  large  side  forces  that 
create  large  rblling  moments  on  the  tow  vehicle.  Although  the 
sonar  does  not  scan  during  a  turn,  it  is  still  necessary  to  keep  the 
towfish  safely  off  the  sea  floor  and  to  minimize  the  above  effects 
while  achieving  steady  state  conditions  on  a  reciprocal  heading 
as  quickly  as  possible.  Large  turning  diameters  reduce  drone  side 
force  but  take  longer.  Natural  ‘overshot’  U-turns  are  shown  to  be 
a  good  compromise  between  these  competing  demands. 

2.  DYNAMICS  MODELLING 

2.1  Drone  Dynamics  (DSSP) 

DSSP  performs  a  time  domain  maneuvering  simulation  by  numer¬ 
ically  integrating  the  6  degree-of-freedom  equations  of  motion  for 
an  underwater  vehicle  (Abkowitz,  1969).  The  drone  starts  neu¬ 
trally  buoyant  in  level  trim  with  a  steady  forward  speed.  At  arbi¬ 
trary  time  intervals,  control  surface  deflection,  ballasting,  flood¬ 
ing,  and  propulsion  changes  are  applied.  Special  modifications 
were  made  in  DSSP  to  accomodate  differentially  deflecting  aft 
planes  (for  roll  control),  to  implement  DOLPHIN-like  autopilots 
for  depth,  heading,  pitch  and  roll  control  (Seto,  1998b),  to  pro¬ 
cess  and  apply  cable  tensions  to  the  towpoint  (Seto,  1998a),  and 
to  model  the  DOLPHIN  Mk2  hydrodynamic  characteristics  and 
propulsion  system. 

Active  first  order  control  algorithms  are  used  for  pitch,  depth 
(heave),  heading,  and  roll  control.  The  system  responds  to  errors 
in  these  quantities  and  their  rates  (the  first  time  derivative  of  the 
error).  The  aft  planes  control  pitch  and  roll  and  the  foreplanes 
control  depth.  On  older  DOLPHINs,  differentially  deflecting  fore¬ 
planes  were  used  for  roll  control  (Watt  et  al.,  1997). 

Drone  added  mass  (Watt,  1988)  and  viscous  forces  and  mo¬ 
ments  are  modelled  conventionally  in  the  equations  of  motion 
using  both  linear  and  nonlinear  hydrodynamic  coefficients.  The 
coefficients  are  estimated  with  standard  methods  used  by  DREA 
and  are  drawn  from  both  theoretical  and  empirical  means  (Watt 
et  al.,  1997). 
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The  propulsion  model  simulates  the  time  dependent  propul¬ 
sion  forces  ex  erted  on  the  vehicle  throughout  a  maneuver  in  which 
engine  rpm,  vehicle  velocity  and  direction  may  be  continually 
changing.  Included  in  the  model  are: 

1)  thrust/power  as  a  function  of  vehicle  speed,  crossflow  angle. 

and  commanded  and  current  rpm: 

2)  engine  power  limitations  as  a  function  of  rpm: 

3)  engine  rpm  time  response  when  an  rpm  command  change  is 

made. 

The  rpm  iime  response  prediction  is  based  on  a  standard,  sec¬ 
ond  order  control  system  model  (Watt,  1990)  that  is  also  used  for 
modelling  other  control  functions  in  the  vehicle,  such  as  control 
surface  deflection  dynamics. 

DYNTOCABS  predicted  tow  cable  forces  are  applied  to  the 
drone  using  DSSP!s  capability  to  model  arbitrary,  time  depen¬ 
dent,  external  forces  acting  on  the  drone  at  arbitrary  locations. 
However,  a  time  history  of  these  forces  must  be  fully  specified  at 
the  start  of  the  simulation,  making  it  necessary  to  iterate  between 
the  two  programs. 

Sea  state  and  free  surface  effects  are  currently  not  modelled 
in  DSSP.  These  can  be  simulated  by  using  DSSP’s  external  force 
capability,  bu  t  have  not  yet  been  implemented. 

2.2  Cable  Dynamics  (DYNTOCABS) 

The  three-dimensional,  finite  difference,  unsteady  cable  code 
DYNTOCABS,  developed  by  Kamman  (Kamman  et  al.,  1989, 
1990.  1996)  at  Coasted  Systems  Station.  Panama  City,  simulates 
the  cable/towfish  dynamics.  The  cable  is  modelled  by  a  series  of 
rigid  links  connected  in  a  chain  with  frictionless  spherical  joints 
to  simulate  a  flexible  cable.  The  3-component  inertial,  fluid  and 
gravitational  loads  are  assumed  to  be  uniformly  distributed  over 
each  link.  These  loads  are  halved  and  concentrated  at  connecting 
joints. 

DYNTOCABS  models  the  towfish  as  a  six  degree-of-freedom 
body  with  mass  and  inertia.  An  arbitrary  reference  point  on  the 
fish  is  connected  to  the  last  cable  link  by  a  frictionless  spherical 
joint.  The  DYNTOCABS  towfish  model  is  described  below. 

The  parameters  specified  for  each  cable  link  are  length,  di¬ 
ameter,  weight,  buoyancy,  hydrodynamic  loading  functions  (force 
distributions),  initial  steady  state  orientation  angle  guesses,  and 
time  derivatives  of  these  guesses.  The  Folb-Nelligan  hydrody¬ 
namic  loading  functions  (no  side  force)  for  armoured  tow  cable 
were  used  (Casarella  and  Parsons,  1970).  The  boundary  con¬ 
ditions  required  for  DYNTOCABS  include  an  acceleration  time 
history  specification  for  the  tow  point  (at  the  drone). 

DYNTOCABS  proceeds  by  first  carrying  out  a  steady  state 
analysis  on  the  input  configuration  which  results  in  a  correct 
steady  state  cable  profile  and  towfish  attitude.  It  then  propa¬ 
gates  the  maneuver  with  a  nonlinear  time  domain  analysis.  As  it 
is  a  boundary  condition,  the  towpoint  acceleration  specification  is 
rigidly  adhered  to.  The  output  from  DYNTOCABS  includes  the 
tow  point  tensions  as  a  function  of  time.  The  resolution  in  time  of 
this  data  can  be  as  high  as  necessary.  Generally,  an  output  inter¬ 
val  of  about  0.5  second  is  enough  for  a  10  knot  maneuver  though 
the  DYNTOCABS  integration  is  performed  every  0.01  seconds. 


2.3  Towfish  Dynamics  (DYNTOCABS) 

In  DYNTOCABS,  the  towfish  dynamics  are  modelled  through  a 
set  of  55  linear  and  nonlinear  hydrodynamic  coefficients  which 
are  incorporated  in  equations  of  motion  similar  to  those  solved 
by  DSSP.  Minor  modifications  were  made  to  this  model  to  ac¬ 
comodate  additional  nonlinear  viscous  damping  coefficients  and 
additional  active  control  surfaces  required  by  the  current  towfish. 
The  towfish  hydrodynamic  coefficients  were  estimated  using  the 
same  methods  as  for  the  drone. 

The  current  towfish  has  a  large,  active  depressor  which  takes 
the  fish  to,  and  maintains,  depth.  The  fish  also  has  four  active 
aft  planes  that  share  control  authority  for  roll,  pitch,  and  yaw. 
Roll  control  can  be  implemented  either  by  yawing  with  the  rud¬ 
der  (to  keep  the  fish  aligned  with  the  tow)  or  by  generating  a 
torque  through  differential  plane  deflections.  Both  methods  were 
used,  with  emphasis  on  the  former.  PID  controllers  are  used  in 
real  life  and  implemented  in  these  simulations.  The  towfish  is 
approximately  neutrally  buoyant.  It  is  equipped  with  an  inertial 
navigation  unit  which  provides  the  unsteady  roll,  pitch,  and  yaw 
error  and  error  rate  signals.  A  depth  transducer  gives  depth  error. 

A  classical  linearized  stability  analysis  was  applied  to  the  tow¬ 
fish,  assuming  a  fixed,  constant  velocity  towpoint.  This  showed 
the  fish  to  be  stable  in  yaw  and  roll  but  unstable  in  pitch.  The 
analysis  predicted  that  the  instability  could  be  eliminated  by  mov¬ 
ing  the  towpoint  forward  about  15%  of  the  towfish  length.  How¬ 
ever,  since  the  towfish  simulation  uses  active  pitch  and  depth 
control,  pitch  instability  was  never  a  problem  and  the  tow  point 
was  not  moved  for  the  simulations  carried  out  for  this  paper. 

3.  COMBINED  DRONE/TOW  SYSTEM  MODEL 

The  procedures  described  here  are  for  simulating  a  U-turn  with 
a  simple  ‘square  wave’  rudder  deflection  profile.  The  rudder  is 
deflected  to  a  constant  value  6r  at  t  =  0  and  back  again  to  zero 
at  t  =  t,.  These  procedures  will  need  to  be  modified  for  different 
maneuvers.  _ 

The  iteration  between  DSSP  and  DYNTOCABS  is  coordi¬ 
nated  in  a  Maple*  worksheet.  A  combination  of  Maple  procedures 
and  executable  fortran  programs  are  called  from  this  worksheet, 
with  plotting  of  variables  available  within  the  worksheet.  The 
iteration  proceeds  as  follows: 

1)  An  initial  drone  trajectory  is  established  using  DSSP  and  as¬ 
suming  tow  cable  tensions  are  zero.  This  establishes  the  rud¬ 
der  commands  and  their  timing  (the  ‘command  profile’)  that 
give  a  turn  with  a  nominally  correct  lane  spacing  (assuming 
the  drone  is  trying  to  follow  a  survey  grid). 

2)  A  three  dimensional  inertial  reference  frame  representation  of 
the  drone  tow  point  trajectory  as  a  function  of  time  is  fit¬ 
ted  with  cubic  splines  and  differentiated  twice  to  form  the  tow 
point  acceleration  time  history  profile  input  to  DYNTOCABS. 
Second  order  acceleration  profile  splines  are  the  usual  DYN¬ 
TOCABS  input.  DYNTOCABS  was  modified  to  accept  this 
input  at  arbitrarily  high  densities  in  order  to  capture  even 
small  scale  nuances  in  drone  motion  (Seto,  1998a).  The  splined 
tow  point  trajectory  data  is  typically  recorded  every  half  sec¬ 
ond  for  a  drone  speed  of  10  knots. 

3)  DYNTOCABS  is  run  and  outputs  three  component  tow  point 
tension  at  pre-determined  time  increments  (about  every  half 
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second)  along  the  vehicle  trajectory.  This  output  is  splined 
and  transformed  to  the  drone  body  axes,  forming  the  input  to 
the  next  iteration  with  DSSP.  Before  proceeding  with  the  next 
DSSP  iteration,  however,  the  rudder  command  profile  must  be 
adjusted  to  deal  with  the  changed  cable  tensions. 

4)  For  the  first  iteration  only,  the  rudder  deflection  needed  to 
offset  the  tow  is  approximated  from  linearized  versions  of  the 
lateral  force  and  yawing  moment  equations  of  motion: 

r'\Yj (m' x'G  -  N'r)  +  NjW  ~  Q]  +  ~ 

Y’N'6r  -  NiYJr 


towpoint  moving  at  10  knots,  and  a  towfish  maintaining  100  m 
depth,  the  comparison  yields: 

2D  Steady  Current 


Quantity 

State 

Model 

towfish  depressive  force  (N) 

7330 

7490 

DOLPHIN  axial  tension  (N) 

8940 

8960 

DOLPHIN  normal  tension  (N) 

1340 

1250 

towfish/drone  axial  separation  (m) 

269 

270 

This  is  good  agreement. 

5.  RESULTS 


Here,  the  average  lateral  body  force  cable  tension,  Ty.  is  an  ex¬ 
ternal  force  added  to  the  equations  of  motion  at  the  tow  point 
xt .  The  desired  yaw  rate  r'  =  rljU  is  that  from  the  initial 
DSSP  calculation  in  step  1.  The  remaining  terms  are  hydrody¬ 
namic  coefficients.  For  subsequent  iterations,  a  quasi-Newton 
iteration  scheme  is  used  to  iterate  the  rudder  deflection  to¬ 
wards  an  angle  that  gives  the  desired  yaw  rate,  r0: 


6r.i  —  + 


(5r,i-l  —  ^r,i‘— 2  )(rp  —  r,--i) 


r,_i  —  r,_  2 


5)  The  time,  t,.  at  which  the  rudder  is  re-zeroed  to  exit  from  the 
turn  is  chosen  so  the  final  vehicle  heading  is  ipf  =  180°.  This 
proceeds  as  follows:  fj,i  =  fa,o  +  At,  t,s 2  —  <»,  1  +  At,  where 
At  is  nominally  10  seconds  (this  has  little  effect  on  the  final 
result)  and,  for  subsequent  iterations: 


tj  ,i  =  tj.,—1  + 


(*«,.-!  -  f,,,~2)(180  -  1) 

1  -  VV,.- 2) 


6)  This  double  iteration  in  6r  and  ta  only  converges  if  the 
two  variables  are  iterated  separately.  Thus,  the  first  DSSP/ 
DYNTOCABS  6r  iterations  are  carried  through  to  a  con¬ 
verged  trajectory  before  the  first  t,  iteration  takes  place.  The 
iterations  are  redone  with  1  (but  are  faster  because  the  tra¬ 
jectory  does  not  change  prior  to  t  =  t, )  and  redone  a  second 
and  third  time  with  t,t 2  and  ta, 3.  The  t,: 3  iteration  results  in 
an  accurate  ipf  =  180°  heading  because,  with  the  tow  system 
turn  rate  at  t  =  ta  almost  steady  state,  ip f  is  a.  linear  function 
of  ta  as  assumed  in  the  above  formula  for  t<lt . 

About  8  DSSP/DYNTOCABS  iteration  cycles  are  required  with 
ta. 0  and  5  with  each  of  ttll  through  ta> 3.  The  total  iteration  takes 
about  1.5  hours  on  a  Pentium  (166  MHz)  personal  computer. 
Trajectory  convergence  is  assumed  when  the  largest  difference  in 
the  trajectories  from  two  consecutive  iterations  is  about  0.001  m. 


4.  MODEL  VALIDATION 

No  data  exists  for  this  particular  cable  and  towfish  system  yet. 
However,  DOLPHIN/towfish  sea  trials  are  planned  to  acquire 
data  for  validating  the  simulator. 

The  simulation  can  be  compared  with  an  ISER  in-house, 
steady-state,  two-dimensional  empirical/ analytical  model  of  a 
towed  system  (that  does  not  account  for  drone/towfish  interac¬ 
tions).  This  program  has  been  tuned  and  calibrated  over  the 
years  and  is  a  reliable  predictor  of  steady  state  straight  and  level 
flight  tow  system  characteristics.  For  300  m  of  cable  scope,  the 


The  interaction  simulation  was  applied  to  the  DOLPHIN/cable/ 
towfish  system  for  U-turn  maneuvers,  as  shown  in  Figures  2 
through  6.  Results  are  preliminary  but  do  show  how  much  in¬ 
teraction  is  present  and  highlight  the  considerations  involved  in 
towing  a  large  fish. 

Figure  2  shows  the  predictions  for  a  non-interacting  (ie,  0th 
iteration)  drone  and  active  towfish  going  through  a  400  m  diam¬ 
eter  U-turn  at  10  knots  with  300  m  of  cable  and  a  100  m  deep 
towfish.  Towfish  depth  keeping  is  good,  pitch  and  roll  angles  are 
within  3  degrees  of  zero,  and  the  tow-off  angle  at  the  drone  (the 
angle  the  tow  cable  is  pulled  away  from  the  drone:s  vertical  plane 
of  symmetry,  as  seen  in  planform)  is  as  large  as  20  degrees.  Inter¬ 
estingly,  though  the  tow  point  turn  is  complete  in  about  2  min¬ 
utes,  cable  tensions  do  not  reach  steady  state  values  until  another 
minute  has  elapsed,  presumably  because  the  cable  dynamic  time 
constants  are  long  (the  cable  is  still  settling  to  a  new  equilibrium). 

Figure  3  shows  the  converged,  fully  interactive  maneuver  that 
results  from  the  Figure  2  initial  conditions.  The  rudder  angle  re¬ 
quired  to  produce  Figure  3  was  a  factor  of  two  larger  than  that 
used  in  Figure  2.  The  overshoot  at  the  end  of  the  Figure  3  tra¬ 
jectory  results  from  the  tow-off  force  continuing  to  pull  the  drone 
inwards  even  though  the  rudder  has  been  zeroed.  Figure  4  shows 
the  effect  of  the  ta  iterations  on  final  heading. 

It  is  expected  that  the  overshoot  could  be  eliminated  in  this 
trajectory  with  an  appropriately  chosen  continuously  varying  rud¬ 
der  command  profile.  However,  the  overshoot  characteristic  is 
desirable.  It  gets  the  cable  behind  the  drone  more  quickly  than 
the  simple  U-turn,  as  demonstrated  below. 

Other  convergence  schemes  are  possible  but  have  not  been 
pursued  at  this  time.  Presumably  one  can  be  found  which  gives 
identical  lane  spacing  for  the  0th  and  final  iterations.  The  main 
objective  here  was  to  generate  full  interaction  solutions  and  ex¬ 
amine  the  magnitude  of  the  effect.  For  this,  comparisons  between 
0th  iteration  and  converged  trajectories  should  be  made  for  runs 
with  the  same  lane  spacing.  This  is  done  in  Figures  5  and  6. 

A  relevant  operational  survey  issue  is  how  long  it  takes  the 
towfish  to  complete  its  turn  and  begin  scanning  another  swath. 
Figure  5  compares  the  towfish  trajectory  from  the  interactive 
model  (Figure  3a)  to  that  for  a  0th  iteration  U-turn  similar  to 
Figure  2a  but  with  the  same  lane  spacing  as  in  Figure  3a.  The 
iterated  model  gets  the  towfish  on  track  faster  because  it  incor¬ 
porates  overshoot. 

A  comparison  of  0th  iteration  and  fully  interactive  cable  ten¬ 
sions  at  DOLPHIN!s  towpoint  show  small  changes  in  form  and 
magnitude  (Figure  6a,c,e).  These  mainly  reflect  the  difference  in 
time  at  which  the  rudder  is  zeroed.  Some  difference  in  the  lat¬ 
eral  tension  (Figure  6c)  occurs  as  the  interactive  DOLPHIN  goes 
through  the  turn. 
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(a)  DOLPHIN  and  towfish  trajectory 


X  (m) 


(c)  DOLPHIN  cable  tow-off  (-)  and  rudder  angle  ( — ) 


(e)  DOLPHIN  cable  tensions 


Fig.  3  Converged  interactive  simulation 

further  validation  is  required  and  full  scale  trials  currently  under¬ 
way  will  allow  this. 

A  preliminary  examination  of  the  interaction  dynamics  during 
U-turn  trajectories  using  a  simple  rudder  command  profile  has 
been  presented.  These  interactions  have  a  significant  effect  on 
the  tow  system  trajectory  and  its  evolution  in  time.  Modelling 
them  has  led  to  a  natural  way  of  turning  incorporating  overshoot 
which  brings  the  towfish  back  on  track  more  quickly  than  occurs 
with  a  simple  U-turn. 


(b)  towfish  depth 


(d)  towfish  pitch  (-)  and  roll  ( — )  angles 


(f)  DOLPHIN  depth  (-),  pitch  (-.),  and  roll  (-  -) 


ith  rudder  re-zeroed  at  ts  = 168  seconds. 

This  work  has  shown  that  interaction  effects  for  DOLPHIN 
and  its  towfish  are  significant:  they  will  likely  be  larger  with 
longer,  deeper  tows.  Since  the  current  iterative  approach  to  solv¬ 
ing  for  these  effects  is  computationally  intensive,  and  requires  cus¬ 
tomized  convergence  algorithms  that  are  trajectory  dependent,  a 
fully  integrated  drone/cable/towfish  dynamic  simulation  elimi¬ 
nating  iteration  is  desirable.  This  is  particularly  true  when  way- 
point  line-following  modes  and  complex  trajectories  are  consid¬ 
ered. 
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Fig.  4  Converged  interactive  trajectories  for  the  0th  through 
3rd  t,  iterations  employed  to  get  the  Figure  3  results. 
The  -F  symbol  on  a  trajectory  indicates  where  the  rudder 
was  re-zeroed. 


Fig.  5  A  comparison  of  357  m  wide  towfish  U-turns. 

Nevertheless,  the  current  approach  can  be  successful  for  simu¬ 
lating  the  simple  trajectories  common  in  route  surveying.  Indeed, 
work  is  underway  to  use  the  current  method  to  assess  the  tow- 
fish  response  to  drone  surge  induced  by  heavy  sea  states  during 
straight  line  surveys. 
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Fig. 


(a)  axial  DOLPHIN  tension 


(c)  lateral  DOLPHIN  tension 


(e)  normal  DOLPHIN  tension 


Cable  tension  levels  at  DOLPHIN’S  towpoint  and  the  towfish 
and  0th  iteration  ( — )  tension  components. 


(b)  DOLPHIN  cable  tow-off  angle 


(d)  towfish  axial  tension 


(f)  normal  towfish  tension 


the  357  m  turns  shown  in  Figure  5:  converged  interactive  ( — ) 
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